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We show how recent data from observations of the cosmic microwave background may suggest the 
presence of additional radiation density which appeared after big bang nucleosynthesis. We propose 
a general scheme by which this radiation could be produced from the decay of non-relativistic matter, 
we place constraints on the properties of such matter, and we give specific examples of scenarios in 
which this general scheme may be realized. 



I. INTRODUCTION 

The era that follows big bang nucleosynthesis and pre- 
cedes matter domination is usually assumed to be domi- 
nated by radiation [T] . The relativistic degrees of freedom 
during that era are also posited, by continuity, to be the 
same as those that prevailed during big bang nucleosyn- 
thesis. However recent measurements by WMAP [5] and 
ACT [3] show a central value for the number of relativis- 
tic degrees of freedom during this era to be higher than 
the one inferred from big bang nucleosynthesis. With 
current data, the difference is not yet statistically signifi- 
cant. In this paper we will take as a working assumption 
that the central value reflects an accurate measurement 
of the relativistic energy density during that epoch. 

There are stringent bounds from big bang nucleosyn- 
thesis, particle physics, and the cosmic microwave back- 
ground on the number and type of light degrees of free- 
dom present in the early universe. These bounds con- 
strain the possible degrees of freedom that might be re- 
sponsible for the new data. Attempts have been made 
to accommodate the new data by introducing relativistic 
sterile neutrinos [4] , or by allowing for neutrino asymme- 
try [5J|S]. We will develop in this paper an alternative ex- 
planation by postulating the existence of a degree of free- 
dom that behaves as non-relativistic matter that is neu- 
tral under the visible gauge group SU(3)cxSU(2)xU(l). 
This degree of freedom can be a non-relativistic particle 
or possibly a coherent oscillating (pseudo) scalar field, it 
contributes a fraction of the energy density during big 
bang nucleosynthesis, and decays into dark radiation af- 
ter light nuclei have been formed. The resulting radia- 
tion adds a relativistic contribution to the energy density 
during the era of interest. 

If future data indeed produces a statistically significant 
difference between the effective number of neutrinos dur- 
ing big bang nucleosynthesis and the number obtained 
from cosmic microwave background observations, then 
this will pose a significant challenge to the sterile neu- 
trino and neutrino asymmetry proposals in their present 
form. 



The paper is organized as follows: Section [TT] will de- 
scribe the relevant observational data which constrain the 
radiation content of the universe. In Section |1Tl| the fo- 
cus is on the implications for the cosmic microwave back- 
ground from the presence of an invisible non-relativistic 
particle decaying after big bang nucleosynthesis into dark 
radiation. We will discuss various regions of parame- 
ters, constrain the properties of the decay products, and 
comment about possible shortcomings of the sterile neu- 
trino or neutrino asymmetry proposal. In Section [rv] we 
present two specific examples where this general proposal 
can be realized. We close with conclusions in Section W\ 



II. OBSERVATIONAL CONSTRAINTS 

It has long been known that primordial light element 
abundances can be used to constrain the energy density 
of relativistic degrees of freedom during the period of big 
bang nucleosynthesis [7]. More recently, observations of 
temperature fluctuations in the cosmic microwave back- 
ground have begun to provide an independent constraint 
on the radiation energy density at a much later period 
in the history of the universe. Specifically, for a fixed 
matter density, increasing the radiation density will re- 
sult in a later time for matter-radiation equality, and this 
has the effect enhancing the early-time integrated Sachs- 
Wolfe effect. Additionally, relativistic particles which are 
not coupled to the photon-baryon fluid are capable of 
free-streaming out of gravitational potential wells faster 
than the sound speed which provides a phase shift and 
damping of the acoustic peaks [8] . 

The bounds on radiation energy density are usually 
given in terms of a parameter N e f f describing the effec- 
tive number of neutrino species defined by 
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The number of light neutrino species is constrained by 
measurements of the decay width of the Z boson [9] to 
be 2.984 ± 0.008. The standard model prediction with 
three flavors of neutrinos and allowing for residual heat- 
ing of the neutrino fluid due to electron-positron anni- 
hilation and other subtle effects [10] is N e ff — 3.046. 
Observations of a value for N e ff which is different from 
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this prediction imply a departure from the standard sce- 
nario of a nearly thermal distribution of three flavors of 
standard model neutrinos. 

Observations of primordial 4 He abundance provide the 
best constraint on N e ff during big bang nucleosynthe- 
sis. These observations are consistent with the standard 
model predictions giving a constraint (at 68% CL) [TTj 
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Some authors have recently reported a higher value for 
the primordial helium abundance |121 113] . The authors of 
[T2"] argue that this may imply a larger value for N B ^ BN , 
while the authors of [13] argue that the uncertainties are 
too large to suggest any tension with the standard big 
bang nucleosynthesis scenario. This is an important is- 
sue which warrants further research, however, for the re- 
mainder of this paper, we will be using the constraint 
quoted above which seems to be consistent with earlier 
studies [T4"] . 

The current constraints from observations of the cos- 
mic microwave background are significantly weaker than 
those from big bang nucleosynthesis. The WMAP satel- 
lite [5] derives its constraints on N e f t primarily from the 
first and third acoustic peaks, while the ground-based At- 
acama Cosmology Telescope [3] uses observations of the 
third through the seventh peaks. These complementary 
measurements which span a broad range of scales give 
some confidence that the results of the measurements are 
not due to a feature which affects only a portion of angu- 
lar power spectrum. The current constraints are (at 68% 
CL) 
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4.56 ±0.75. (4) 



The constraints listed here assumed ACDM cosmology 
with massless neutrinos. Allowing for non-zero neutrino 
masses or varying dark energy equation of state affects 
the allowed parameter range [4] [15] . 

These limits are consistent with N B f B at the 2-cr level, 
but there is a hint of tension between the central values 
of these measurements. The Planck satellite will signif- 
icantly increase the precision of cosmic microwave back- 
ground observations [16] giving 5N^ B ~ 0.26 or better 
[17j , and it will also help to break degeneracies with non- 
zero neutrino masses and varying equation of state for 
dark energy [15] . If the central values from the current 
measurements reflect an accurate measurement of N e ff, 
Planck could reveal a disagreement between N B j BN and 

N°/f B at the level of 4- to 5-cr. 



III. DECAYING MATTER 

If Planck reveals a significant discrepancy between 
N B * BN and N^Jf B , it is important to ask how these re- 
sults may be reconciled. We suggest that the decay of 



non-relativistic matter into radiation would provide the 
necessary increase in N e f f to resolve the conflict. The en- 
ergy density in the relativistic decay products, and hence 
the contribution to N^ IB , will be determined by the en- 
ergy density in the non-relativistic matter at some earlier 
time along with the lifetime of the matter. For matter 
with an energy density px at t — 10~ 4 s, and a lifetime 
r, we can calculate the increase in N e ff after the decay 
to be 
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where po is the photon energy density at t = 10 4 s 
adjusted to include the effect of heating due to electron 
positron annihilation, and is given by po = or (1.334 x 
10 12 K) 4 = 1.149x 10~ 4 GeV 4 Q]. In deriving {Hf, we have 
used the fact that the energy density of non-relativistic 
matter decreases as a~ 3 and that of radiation as a -4 , 
where a is the Robertson- Walker scale factor which scales 
as t 1 / 2 during the radiation dominated era. We have also 
assumed a very rapid decay into radiation at the time 

t = T. 

It is important that this additional matter give at most 
a sub-dominant contribution to the total energy density 
during the period of big bang nucleosynthesis so as not 
to alter the predictions of primordial light clement abun- 
dances. We thus require that N e ff < 3.2 for t < 20 
minutes. For a more detailed study of the constraints on 
decaying particles during big bang nucleosynthesis see 
[15] [T9] . It is also necessary that the decay into radia- 
tion takes place before the observable modes of the cos- 
mic microwave background have reentered the horizon. 
This requirement constrains t < 1650 years [20]. Putting 
this together, we can find the allowed range of parame- 
ters which satisfies the constraints from both N B f BN and 
N^ff B ■ We summarize our results in Figure [TJ 

It seems necessary that the relativistic decay products 
be neutral under the standard model gauge group. Light 
particles carrying strong or weak charge are very well con- 
strained by collider experiments [5]. A decay into pho- 
tons would run the risk of destroying the nuclei produced 
during big bang nucleosynthesis, and thus altering the 
primordial light element abundances. Furthermore, these 
photons would equilibrate in the photon-baryon plasma, 
increasing the photon temperature relative to the neu- 
trino temperature, and therefore decreasing rather than 
increasing N e ff as measured in the cosmic microwave 
background. This means that if a scenario such as the 
one we have described is realized in nature, it is unlikely 
that we would ever be able to observe the effects of the 
new contribution to the energy density, except through 
its gravitational signature. 

A suggestion similar to the one made here appeared 
previously in the literature which was motivated by data 
from large scale structure and Lyman-a forest measure- 
ments [H] . The cosmological implications of that model 
were discussed in detail in 122 . 
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(a)WMAP 7- Year Constraints 
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(b)Planck Projected Constraints 

FIG. 1. These figures illustrate the constraints on energy den- 
sity and lifetime of non-relativistic matter decaying into radia- 
tion. The vertical axis gives the ratio of the energy density with 
po = Ob(1.334 X 10 12 K) 4 = 1.149 X 10" 4 GeV 4 evaluated at 1CT 4 
seconds, and the horizontal axis shows the lifetime reported in sec- 
onds. The blue shaded region near the top of the figure is excluded 
due to the constraint that N e ft < 3.2 during big bang nucleosyn- 
thesis while the shaded region near the right is excluded by re- 
quiring that the decay into radiation occurs before the highest I 
modes observable in the cosmic microwave background reenter the 
horizon. The dark and light orange shaded regions give the 1- and 
2-a constraints, and the black dashed line corresponds to the cur- 
rent observed central value of N C M B =4.34 from WMAP. The 
constraints from WMAP 7-year data are consistent with no extra 
energy density, while the projected constraints from Planck require 
a nonzero energy density and fix the lifetime of the decaying matter. 

We have focused on the case where the energy density 
before the decay is in the form of non-relativistic mat- 
ter, however, this is not the only possibility. We require 
only that there is some component of the energy den- 
sity which redshifts more slowly than radiation so that 
the ratio px/p-y increases with time. This allows px to 
be subdominant during big bang nucleosynthesis while it 
may become important before the modes relevant for the 
cosmic microwave background have reentered the hori- 
zon. Note that scenarios such as the addition of sterile 



neutrinos [3], other decoupled radiation [33], or neutrino 
asymmetry [3 [6] tend to give N® f B f N = N^/f B , and thus 
if future data reveals a statistically significant difference 
between N B j BN and N^ IB , these models would face 
a significant challenge. The requirements on the decay 
products are more restrictive because they must act as 
radiation in order to have the same effect as an increase in 
N e ff. These particles need not be strictly massless, but 
it is important that they are relativistic until well after 
matter-radiation equality so as not to significantly alter 
late time cosmology [23]. For a fluid component with a 
pressure given by px = wpx with w < 1/3, equation ^ 
becomes 

— - G)("f d^r* 2 ^™ 
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IV. EXAMPLES 
A. Non-Relativistic Invisible Particle 

A context in which such non-relativistic particles oc- 
cur naturally are supersymmetric models which exhibit 
a hierarchy of splittings in super- multiplets 25j. This 
hierarchy was exemplified in [25] in the context of gauge 
mediation [26H35] with mediators being charged under 
the standard gauge group as well as an invisible gauge 
group. The hierarchy in masses of the mediator fields 
reflects itself in the hierarchy noted in splittings. Details 
of the construction can be found in the paper by Fis- 
chler and Tangarife-Garcia [S3]- For other models with 
decaying particles see for example [5TJ [35] I31)H3"5] . 

A prime candidate as shown in the aforementioned pa- 
per is an invisible photino. This particle is neutral under 
the visible gauge group, SU(3)c x SU(2) x U(l) and is 
the super-partner of an invisible U(l) gauge field. The 
low mass of this particle is a consequence of the small 
splitting inherited by the super-multiplets in the invisi- 
ble sector. The invisible photino decays into an invisible 
photon and a relativistic gravitino. It's lifetime, cross- 
sections, and mass are related to various parameters in 
the class of models mentioned above. These parameters 
include in particular, the supersymmetry breaking scale, 
the mass of the relevant mediators, and gauge coupling 
constants. The abundance of the invisible photino is fixed 
by the property that it is in equilibrium after reheating 
with other invisible degrees of freedom including invisible 
photons. As noted in |25j the temperature of the invisible 
sector is lower than the visible sector. Below, we exhibit 
the values for quantities related to parameters occurring 
in a specific model belonging to this class and adjusted so 
that they are in agreement with the values obtained from 
Figure [l] An additional requirement is that the gravitino 
mass not affect late time cosmology. 

In particular, we consider three parameters: the SUSY 
breaking scale p, the photino mass m^,, and the annihi- 
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lation cross section (crv). These are given as functions of 
t by the following expressions. 
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x f 



2.098 x 10 15 GeV^ 



where /3 is the ratio between the current radiation energy 
density and the current gravitino energy density 



Pg,o 



and Xf = m^/Tf is the "freeze-out" parameter which is 
commonly taken to be i/ ~ 20. For the specific model 
presented in the Fischler-Tangarife paper, these expres- 
sions can be rewritten in terms of a new set of parameters 
77i2, A* an d Vh by using the relations: 
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The above generic constraints then translate into re- 
strictions on the value of M,m, <fio, where M is the de- 
nominator mass of the dimension 5 operator, m is the 
mass of the (pseudo)scalar field, and 4>o is the value of <f> 
at the onset of oscillations. 

The parameters are given by the following expressions 



m (r) = (2.417 x 10~ 27 GeV)£ 



Is 



M(r) = (1.188 x 10~ 28 GeV)^ 6 ( " ) AiV^/ , (16) 



AN 2 ef} (15) 



where £ defined by 

4>o = £M P , 

can take values in the range 

3 x 10~ 9 < £ < 10~ 7 . 

For the upper bound of £ (£ ~ 10~ 7 and using N e jf ~ 
4.6), the minimum values that can be obtained for 
and M are: 

m ~ 10 _1 eV, M ~ 10 3 GeV . 

For the lower bound, ^ ~ 3 x 10 -9 , we get that the 
maximum values for these masses are: 

~ 10 2 GeV, M ~ 10 18 GeV . 



B. Oscillating (Pseudo) Scalar 



CONCLUSIONS 



Another possible realization of our generic proposal is 
to consider a (pseudo) scalar field whose homogeneous co- 
herent oscillations mimic non-relativistic matter and con- 
tribute a sub-dominant component of the energy density. 
This (pseudo) scalar field decays into dark photons via 
the following dimension 5 operator: 
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The potential for this (pseudo) scalar field is assumed to 
be harmonic in the region of field space where the co- 
herent (pseudo) scalar field is located at the onset of its 
oscillations. 
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In a nutshell, the (pseudo)scalar field starts oscillat- 
ing when the Hubble constant is of order of the mass of 
the (pseudo) scalar field. The lifetime can be obtained 
from the dimension 5 operator mentioned above. It is 
then straightforward to find how the energy density of 
the background at the time of decay, p T depends on the 
energy density of the background at the onset of oscilla- 
tions po- 
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If one assumes the effective number of neutrino species 
before and after big bang nucleosynthesis is indeed ac- 
curately given by the central values of the experimental 
data, a conundrum ensues requiring the existence of a 
new kind of dark radiation which appears after the for- 
mation of light nuclei. With this assumption current re- 
alizations of sterile neutrinos and neutrino asymmetries 
cannot account for the difference between N~/f N and 
N^^ IB , however with the present experimental precision 
these scenarios are consistent within the experimental er- 
rors. We have shown how the decay of non-relativistic 
matter into dark radiation after big bang nucleosynthe- 
sis may reconcile a difference in N B M N and N^Jf B if 
forthcoming data indeed reveals such a difference. Our 
results for the constraints on the energy density and life- 
time of such matter are summarized in Figure [T] We 
described two specific examples in which this general sce- 
nario could be realized: a supersymmetric model with a 
hierarchy of splittings in supcr-multiplcts, and an oscil- 
lating (pseudo) scalar field. 
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